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Properties of low-dimensional spin-Peierls systems are described by using a one dimensional S = ^ 
antiferromagnetic Heisenberg chain linearly coupled to a single phonon mode of wave vector vr 
(whose contribution is expected to be dominant). By exact diagonalizations of small rings with 
up to 24 sites supplemented by a finite size scaling analysis, static and dynamical properties are 
investigated. Numerical evidences are given for a spontaneous discrete symmetry breaking towards 
a spin gapped phase with a frozen lattice dimerization. Special emphasis is put on the comparative 
study of the two inorganic spin-Peierls compounds CuGeOs and NaV205 and the model parameters 
are determined from a fit of the experimental spin gaps. We predict that the spin-phonon coupling is 
2 or 3 times larger in NaV2 05 than in CuGeOa. Inelastic neutron scattering spectra are calculated 
and similar results are found in the single phonon mode approximation and in the model including a 
static dimerization. In particular, the magnon 5 = 1 branch is clearly separated from the continuum 
of triplet excitations by a finite gap. 

PACS numbers: 64.70.Kb, 71.27.+a, 75.10.Jm, 75.40. Mg, 75.50.Ee 
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I. INTRODUCTION 

Recently, a renewed interest for one dimensional (ID) 
spin chains was created by the observation of spin-Peierls 
transitions in the inorganic compounds CuGeOs ||l| and 
Q!'-NaV205 (|-||]. Below some critical temperature Tspi 
the spin-Peierls phase is experimentally inferred from a 
rapid drop of the spin susceptibility. The low tempera- 
ture phase is characterized by the opening of a spin gap 
A (see below) and the dimerization of the lattice along 
the chain direction as confirmed for example by X-rays 
diffraction in CuGeOs § and NaVaOg |] or by Na NMR 
experiments in NaV205 Q. 

In general, these compounds are well described above 
the transition temperature by a ID frustrated antiferro- 
magnetic (AF) Heisenberg chain. The nearest neighbor 
and next-nearest neighbor spin exchange integrals J and 
J' can be determined by a fit of the magnetic suscepti- 
bility X at high temperatures. In fact the position of the 
maximum of the curve and more generally the magnetic 
properties only depend on the frustration ratio a — J' / J . 

Parameters such as J = 160 K and a — 0.36 ||^ or 
J = 150 K and a = 0.24 S have been suggested for 
CuGeOa. Further studies PjlQ] seem to confirm that the 
dimerization is large in this system and we shall take 
a = 0.36 in the rest of the paper. For the novel NaV2 05 
system J = 440 K and a « have been proposed (5| in 
good agreement with Refs. [p|jll|]. 

The zero temperature spin gap A has been deter- 
mined by several means. Inelastic neutron scattering 
(INS) gives a direct measure of it. So far INS has 
been performed on single crystals of CuGeOa and pow- 
der NaV205 samples. Values of A ~ 2.1 meV [F2|Jl3| 



and A = 9.8 meV ^ have been reported, respectively. 
Other more indirect methods like NMR can also pro- 
vide a measure of the magnitude of the spin gap. ^^Cu 
or ^'^Cu NMR have been performed on single crystals of 



CuGeOa 
NaVaOs 



|14| and ^'^Na NMR on aligned polycrystals of 
4|7 The local magnetic susceptibility is propor- 
tional to the NMR Knight shift and the spin gap is esti- 
mated by a fit of the temperature dependence of the local 
susceptibility below the transition temperature. A value 
of A ~ 8.4 meV was found for NaV2 05. A third indepen- 
dent estimation of A can be obtained by a fit of the low 
temperature bulk magnetic susceptibility of single crys- 
tals measured for instance by a SQUID-magnetometer 
technique. Ref. |^ reports a value of A ~ 7.3 meV for 
NaVsOg. 

The dimensionless ratio A/ J is the crucial parameter 
needed in our theoretical analysis. Results for CuGeOa 
are now well established and the value of A/ J — 0.151 
is often used in the literature. Nonetheless ratios such 
as A/J = 0.203 §, A/J = 0.175 |] or A/J = 0.193 § 
can be found for NaV2 05. Since this last estimation was 
obtained from experiments performed on single crystals 
we have thus decided to use it as a reference. In any case, 
the quite small differences between the previous experi- 
mental values are not relevant. 

Theoretically, the spin dynamics of the ID Heisenberg 
chain depends strongly on the frustration parameter a. 
Indeed, for a > ttc — 0.241 a gap appears in the spin 
excitation spectrum jl^,|l^. Therefore, we expect that 
the two previous spin-Peierls compounds will have quite 
different magnetic properties. CuGcOa is dominated by 
intrachain frustration. On the other hand, NaV205 will 
behave, at high temperatures, more closely to an unfrus- 
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trated Heisenberg chain and, at low temperatures, the 
small interchain frustration alone cannot be responsible 
for the opening of a spin gap. The coupling to the lat- 
tice is therefore expected to play a dominant role in the 
transition at least for NaV205. In order to study their in- 
terplay, the frustration and the spin-lattice coupling have 
to be treated on equal footings. This is the purpose of 
this paper. 

It is well known that a ID system shows no phase tran- 
sition at finite temperature because of quantum fluctu- 
ations. Interchain couplings are necessary to obtain a 
finite transition temperature. However, they are thought 
to be small and will be neglected hereafter in the study 
of zero temperature properties. 

So far, there have been various attempts to treat the 
coupling to the lattice by considering a static dimeriza- 
tion S of the exchange integral (so called adiabatic ap- 
proximation or frozen phonon approximation). The value 
of S is determined in order to obtain the experimental 
value of the zero temperature spin gap A. Dimeriza- 
tions such as 6 = 0.014 § and (5 = 0.048 were 
proposed for CuGeOa and for NaV205, respectively, in 
order to reproduce the measured spin gaps (assuming 
A ~ 0.151J and A ~ 0.193J for CuGeOg and NaVzOs, 
respectively). Calculations using this approach have been 
performed in order to make first comparisons with exper- 
iments 



In this paper we use a modification of the previous 
static model to describe the physical properties of one- 
dimensional spin-Peierls compounds below the transition 
temperature. For convenience, the previous ad hoc static 
dimerization discussed above is replaced here by a sin- 
gle dynamical optical phonon mode (Section ||). As far 
as thermodynamic properties are concerned, this model 
should be, in fact, equivalent in the thermodynamic limit 
to a model where the lattice is treated at a mean field 
level [pl| . However, this new approach has some advan- 
tages: (i) it incorporates automatically the elastic en- 
ergy and avoids the lengthy iterative procedure needed 
in a mean-field treatment to converge to the equilibrium 
static lattice dimerization; (ii) it enables to study the 
mechanism of the lattice symmetry breaking and, hence, 
provides a basis for future studies including a macro- 
scopic number of phonon modes (i.e. proportional to the 
system length L) ||2^ in spin-Peierls chains. 

Within this single mode approximation, we truncate 
the Hilbert space of the phonons and show in details that 
this approximation is well controlled (Section |l|). Us- 
ing a finite size scaling analysis (discussed in detail in 
Section IV) the dimerization and the spin gap resulting 
from a spontaneous discrete symmetry breaking of the 
lattice periodicity are calculated (Section 0). Focussing 
primarily on CuGeOs and NaV205 materials, we then 
establish a simple relation between the parameters of the 
model in such a way to enforce the constraint that the 
numerically calculated spin gap is equal to the experi- 
mental gap. The role of the parameters is discussed. In 
the last part (Section VI), we study the spin dynamics. In 



particular, we investigate the role of the lattice dynam- 
ics on the low energy magnon branch and low energy 
structures in the dynamical spin structure factor. Our 
results are compared to the ones obtained in the static 
model IT7l|T9|^,ra. 



II. MODELS 

Our starting point is the ID frustrated AF Heisen- 
berg chain. For practical applications, the previous val- 
ues of J and a will be used : J = 160 K, a = 0.36 
for CuGeOa § and J = UO K, a = for NaVaOs [|. 
In addition, a coupling between spins and dispersionless 
optical phonons (magneto-elastic coupling) is considered. 
For sake of simplicity we assume a linear dependence of 
the exchange integrals on the relative atomic displace- 
ments {u,} |2|,|§, 

H = Jj2iil + Xu,) S,.S,+i + aS,.S,+2) + i?ph, (1) 



where A is the coupling constant. H^^ is the phononic 
Hamiltonian of identical independent quantum oscilla- 
tors, = + {Pi is the conjugate mo- 
mentum associated to the atomic displacement Ui). The 
atomic displacements Ui and their conjugate variables 
can easily be expressed in term of the canonical phonon 
creation and annihilation operators 6^ and b^- Since the 
spin susceptibility diverges (for a < 0.5 |27 2^) at mo- 
mentum k — TT we expect that the coupling to the lattice 
will be dominant at A; = tt which corresponds, in fact, to 
the modulation of the spin-Peierls ground state. There- 
fore, from now on, we shall only keep a single k = tt 
phonon mode p9]. In this case, using. 



(fi^ = K/m and L is the number of sites), the final 
Hamiltonian becomes. 



(2) 



where g 



2mn 



is the dimensionless coupling con- 



stant. Within this approximation can be rewrit- 
ten as ifpjj = [b\b^ + ^) where VL is the energy of a 
phononic quantum. 

Before going further, we can already discuss qualita- 
tively the physics contained in Hamiltonian (^). Indeed, 
we expect in the thermodynamic limit a discrete symme- 
try breaking corresponding to a doubling of the unit-cell. 
This can be described very simply at the MF level. By 
assuming a dimerization 5 given by the order parameter 



2 



{b^ + h\)MF and omitting a constant part, the 



g 

Hamiltonian takes the form, 



Hmf = J^((l + (5(-ir)5'..S',+i +aS,.5,+2) 



(3) 



case of CuGeOa and a couphng constant g = 0.5 is used) 
for chains of length L — 12 and 20 plotted as a function of 
-^max- Typically, the mean occupation number is smaller 
than 3 as in Figs. || and ^ (and in fact even smaller than 
^ 0.5 for more realistic parameters) and the energy has 
converged for A^max ^ 30. Fig. |^ proves that the trun- 



where the last term is the elastic energy loss. This is ex- 
actly the well known model describing a static dimeriza- 
tion below the transition temperature in spin-Peierls sys- 
tems . Interestingly enough, a similar effective model 
has also been used to describe conjugated hydrocarbons 
with bond alternation such as polyacetylene |^^. In this 
new form, the breaking of the lattice periodicity is ex- 
plicit. As a consequence the ground state becomes doubly 
degenerate (the order parameter d can take a positive or a 
negative value) and a spin gap appears. The spin-Peierls 
ground state is characterized by a ■ ■ A— B — A — B ■ ■ 
pattern with a succession of strong singlet A bonds and 
weak singlet B bonds (so called Valence Bond or dimer 
state). Note that 5 in model (^) is a variational param- 
eter to be determined in order to minimize the ground 
state energy by an iterative procedure. In contrast, the 
dimerization in Hamiltonian (^ arises from a dynamical 
symmetry breaking. However, it is interesting to notice 
that models (||) and (^ should be in fact equivalent ]3l| ] 
in the thermodynamic limit, at least as far as their ther- 
modynamic properties are concerned |21 3^ . 

Static and dynamical quantities are given by exact di- 
agonalizations of small chains. Using a finite size scaling 
analysis, results in the thermodynamic limit are deduced. 
The parameters S on one hand and g and fi/J on the 
other hand are determined from a fit to the experimental 
spin gap. 



III. TRUNCATION PROCEDURE 

Let us now deal first with the numerical treatment of 
(^. The total Hilbert space can be written as the ten- 
sorial product of the space of the spin configurations (to 
which the symmetry group of the problem is applied) 
times the phononic space. However, strictly speaking, 
the Hilbert space associated to the phonons is infinite 
even for a chain of finite length. Indeed, the natural basis 
{\n)} is simply defined by the unlimited occupation num- 
ber n of the fc = TT phonon mode, |n) = l/\/nI(f't)"|0)- 
Such a difficulty can nevertheless be easily handled in 
an exact diagonalization treatment . The solution is 
to truncate the phononic space so that the occupation 
number is smaller than a fixed A^max which has to be 
chosen in an appropriate way. Clearly, if iVmax is, let us 
say, an order of magnitude larger than the exact mean 
occupation number {b^j^b^) the truncation procedure will 
not affect the accuracy of the results which can then be 
considered as basically exact. This can be seen in Fig. |l| 
showing the energy per site of the ground state in the 
spin and 1 sectors (the value of a corresponds to the 
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FIG. 1. Convergence of the energy per site of the lowest 
singlet (o, •) and triplet (□, ■) states in units of J as a func- 
tion of the maximum number of phonons A^max. Parameters 
are a — 0.36, g = 0.5 and Q = 0.3J. Open (filled) symbols 
correspond to L = 12 (L = 20) sites. 

cation procedure is very well controlled even for rather 
(unphysically) large coupling constants like g = 0.5. The 
results reported in the rest of this paper are then ob- 
tained with a sufficiently large value of A^max and the 
preliminary studies of the convergence of the results with 
increasing A^max, although not mentioned each time, have 
been performed for each choice of the parameters of the 
model. 

It is interesting to study the dependence of the mean 
occupation number on the three parameters of the prob- 
lem (length of the chain L, coupling constant g and fre- 
quency of the phonons fl) since, first, this number di- 
rectly determines the practical value of A^max to be chosen 
and, secondly, it provides some physical understanding. 
Fig. |l| already suggests that the mean occupation num- 
ber increases with the length of the chains. To investigate 
this effect in more details, the mean occupation number 
(*ol4&^l*o) in the ground state I'^o) is plotted in Fig. || 
as a function of L. Clearly, (5'o|6jr^7r|^'o) (as well as the 
value required for A^max) grows linearly with the chain 
length L. In fact, this effect is directly connected to the 
breaking of the lattice symmetry as can be seen very 
easily from a very simplified version of Hamiltonian (|^). 
In a symmetry broken state, an effective (approximate) 
phononic Hamiltonian Hph can be constructed by tak- 
ing MF values for the spin operators. Assuming that 
J2ii~'^)''{^i-^i+i)MF (dimer order parameter) varies lin- 
early with L one then gets Hph — Ag^/L{bTr+bl)+^ bl^b^^ 
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(A is an undetermined constant). In this approximation, 
(fet ^tt) = A'^g'^L/il'^ grows linearly with the length of the 
chain. In addition, this simple argument also suggests 
that the occupation number of the tt mode scales like 
the square of the dimensionless coupling g and like the 
inverse square of the phonon frequency. These intuitive 
behaviors are indeed well followed as can be seen in Fig. ^ 
in a large range of parameters. 



0.20 



A 

+ 

V 



0.15 



0.10 



0.05 - 



0.00 















O' 



5 10 15 20 25 

L 



FIG. 2. Dependence of the mean occupation number on 
the length of the chain L for g — 0.109 and Q — 0.3 J. 



formed to obtain extrapolations to the thermodynamic 
limit. 



IV. FINITE SIZE SCALING ANALYSIS 

Firstly, we focus on the size dependence of the en- 
ergy per site of the singlet ground state and of the low- 
est triplet state which are expected to converge to the 
same value in the thermodynamic limit. Typically, we 
use chains of length L = 8, 12, 16, 20 and 24 sites. 
Data are shown in Fig. || for a = 0.15, g = 0.45 and 
= 0.3J. The ground state energy per site varies 
roughly like This behavior is predicted for gapless 

ID chains obeying conformal invariance but seems 
to be still valid here in spite of the presence of a spin gap 
(see later). This already suggests that, for such param- 
eters, the system sizes are still comparable to the spin 
correlation length but not much larger. The behavior of 
the triplet energy is more involved. An approximate 
dependence is expected (giving a square root singularity 
in the 1/i^ units of Fig. ^ if there is a finite spin gap 
A (defined by the L ^ oo extrapolation of the difference 
A(L) = Eo{S = l,L)-Eo{S = 0,L) of the toia/ energies 
of the lowest states of the singlet and triplet spin sectors). 
Such a behavior seems indeed to be observed in Fig. ^. 
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FIG. 3. Mean occupation number calculated on a L = 12 
site chain versus for Q,/J = 0.3 (a) and versus (J/Q)^ for 
g = 0.109 (b). 

One also observes in Fig. |l| that the singlet ground state 
energy is almost converged for L = 20 (the values of the 
energies for L = 12 and L = 20 at large A^max are indis- 
tinguishable) while finite size effects are still large for the 
triplet energy because of the existence of a continuum 
of states above the first triplet excitation. In the next 
Section, we investigate carefully the convergence of var- 
ious physical quantities with respect to the system size. 
We show that an accurate finite size analysis can be per- 
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FIG. 4. Convergence of the energy per site of the ground 
state of spin (o) and 1 (□) in units of J as a function of the 
inverse of the square length of the chain 1/L^ for a = 0.15, 
g = 0.45 and D. = 0.3J. 

Let us now examine in details the behavior of the spin 
gap A(L) versus L to extract values in the thermody- 
namic limit. Requiring that the extrapolated ratio A/ J 
of the model (||) is equal to the observed experimental 
value will lead to some constraint on the model param- 
eters n and g. Our procedure can be summarised in 
three steps; (i) a controlled truncation procedure of the 
phononic Hilbert space for a large set of parameters g, 
f2 and system sizes L, (ii) a finite size scaling analysis 
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in order to accurately determine the spin gap as a func- 
tion of g and 17; (iii) a determination of the relation to 
be followed by the parameters g and il in order that the 
calculated ratio A/ J equals the actual experimental ratio 
(see Section 0) . 

We first consider the scaling behavior of the spin gap. 
We have found that it scales accurately according to the 
law, f^lToll 



0.6 



A(i) = A + yexp(-^), 



(4) 



where Lq is a typical length scale. In general Lq is of 
the order of the magnetic correlation length character- 
izing the decay of the equal time spin-spin correlation 
in real space. As seen later, values of Lq are typically 
20 lattice units (l.u.) for parameters corresponding to 
CuGeOa and 30 l.u. for NaV205. Therefore, with chains 
lengths up to 24 sites, finite size effects are still impor- 
tant and an accurate extrapolation is necessary. This 
scaling is illustrated for a — 0.15 < etc, g — 0.22, 
= 0.3J (o), for a = 0.36 > (CuGeOa-like case), 
g = 0.089, ft = 0.3J (□) and for a = (NaVsOg-like 
case), g = 0.40, n = 0.5J (0) in Fig. |( a). A spin gap 
opens for all a if g > 0. This is similar to the mean- 
field treatment where the order parameter 5 7^ leads to 
the symmetry breaking and thus to the opening of a spin 
gap. 

In Fig. ^(b) we compare, in the case of NaV205 (i.e. 
a = 0), the scaling of the spin gaps calculated using the 
dynamical model (^) with g = 0.275, H. = 0.3J (o) on 
one hand and the static model (^) with S = 0.05 (□) on 
the other hand |36| . These values of the parameters have 
been chosen in order to obtain the same extrapolated spin 
gap. Although the spin gaps are equal, the two models 
exhibit slightly different scaling behaviors (Lq — 30 for 
the dynamical model and Lq ~ 18 for the static one pTf). 

At this stage, it is interesting to better understand how 
in the the dynamical model (0) the opening of the spin 
gap is connected to the discrete symmetry breaking (as 
can be seen e.g. in X-rays scattering). The first signature 
of this phenomenon is the degeneracy of the ground state 
which is expected in the thermodynamic limit. We have 
therefore studied the behavior with system size of the en- 
ergies Ep{S = 0), p = 0, 1, 2, of the three lowest singlet 
states. The energy differences Ei{S = 0) — Eo{S — 0) 
(circles) and Ei{S = 0) — Et^i^S = 0) (squares) are plot- 
ted in Fig. ^, in the case = 0.3J, as a function of the 
inverse length of the chain 1/L for a = 0.36 (open sym- 
bols) and for a = (filled symbols). The values of the 
coupling g are chosen here in such a way to reproduce 
the experimental spin gaps of the CuGeOa (open sym- 
bols) and NaV205 (filled symbols) materials (see Section 
V). The results show very convincingly that the singlet 
ground state is indeed two- fold degenerate in the thermo- 
dynamic limit while a finite gap for singlet excitations 
appears above It is important to notice that the 

quantum numbers associated to the translation symme- 
try are different for the two lowest singlet states which 
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FIG. 5. (a) Spin gap A in units of J as a function of the 
inverse of the length of the chain 1/L for a = 0.15, g = 0.22, 
Q. = 0.3J (o), a = 0.36, g = 0.089, D. = 0.3J (□) and a = 0, 
g = 0.40, n = 0.5 J (0). (b) Comparison between the behav- 
iors A/J vs l/L obtained within the dynamical model ^ for 
a — 0, g — 0.275, and Q — 0.3 J (o) and within the static 
model im ior a = 0,5 = 0.05 (□). 



correspond to momenta k = and k = tt. Hence, mixing 
of these two states leads to a doubling of the unit cell. 

The lattice dimerization can be quantitatively mea- 
sured by the order parameter S* — -^{(b^ + 6^ )^)^/^ 

(the expectation value {b^ + b^) vanishes because small 
tunnelling between the two degenerate dimer states al- 
ways exists in a finite chain). S* as a function of the 
inverse length of the chain is plotted in Fig. |^ for various 
pairs of parameters (f2, g) (see caption) chosen in such 
a way that the spin gap is constant (in fact adjusted to 
the actual spin gap of CuGeOs as described in Section 
V) . Extrapolated values of the dimerization 6* for differ- 
ent phonon frequencies are in fact quite close, at least 
in the range 0.1 < < 0.5. The dimerization S* seems 
then to be only determined by the magnitude of the spin 
gap. The fact that S* , at fixed extrapolated spin gap, 
is independent of the frequency Q is consistent with the 
proof by Brandt and Leschke that the thermody- 
namic properties of the dynamical model (^ and of the 
static model (||) are identical. However, it is interesting 
to notice that the value obtained here (~ 0.022) is sig- 
nificantly larger than the value {^-^ 0.014) needed in the 
MF approximation to produce the same gap. The differ- 



5 



0.4 

O 

LU 

I 

LU 

^ 0.2 

o 

LU 

V 0.1 
LU 



,o''' 



0.0 
0.00 



,0' 



0.05 0.10 0.15 

1/L 



FIG. 6. Energy differences {Ei{S = 0) - £'o(5' = 0))/J 
(o,.) and (£2(5' = 0) - £0(5' = 0))/J (□,■) as a function 
of 1/L. Open and filled symbols correspond to a = 0.36, 
g = 0.109 and Q = 0.3 J and to a = 0, y = 0.270 and 
Q = 0.3J, respectively. 



ence between these two values can be simply attributed 
to the zero point motion of the harmonic mode which is 
included only in (0). 
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FIG. 7. Order parameter S* as a function of the inverse of 
the length of the chain 1/L for Q — O.IJ and g — 0.062 (o), 
D. = 0.3J and g = 0.109 (□) and !^ = 0.5J and g = 0.141 (0) 
(see text regarding the choice of parameters). 



V. COMPARISON WITH EXPERIMENT 

The systematic finite size scaling described above has 
been performed for a large set of parameters Q/ J and g. 
For simplicity, let us first assume fl/ J — 0.3. The behav- 
ior of A(ri/ J = 0.3, g)/J versus g is plotted in Fig. ^ for 
a large frustration a = 0.36 corresponding to the case of 



CuGeOa (o) and for a non frustrated chain correspond- 
ing to the case of NaV205 (□). Quite generally, the spin 
gap grows with the coupling constant g as expected. In- 
deed, a larger coupling to the lattice produces a larger 
dimerization and then, indirectly, a larger spin gap. 

The actual physical value of the ratio fl/J is, to the 
best of our knowledge, difficult to obtain from exper- 
iment. Therefore, we shall not here restrict to any 
specific value of f2/J but rather consider a wide range 
0.1 < < 0.5. However, for each value of SI, the di- 
mensionless coupling constant g{fl) can be determined by 
enforcing that the extrapolated spin gap ratio A(f2, g)/ J 
equals the experimentally observed gap. The procedure 
is shown in Fig. || for fl ~ 0.3 J and a — (NaV205) and 
a = 0.36 (CuGeOa). The small horizontal marks corre- 
spond to the actual experimental gaps, i.e. A/J ~ 0.151 
and A/ J ~ 0.193 for CuGeOs and NaV205, respectively. 
We then obtain g{n = 0.3) ~ 0.109 for CuGeOs and 
g{n = 0.3) = 0.270 for NaVsOs- The same method was 
performed for two other values of the frequency, SI = O.IJ 
and = 0.5 J. A relation is then obtained between SI and 
g for the two values of the frustration parameter a — 
and a = 0.36. This is illustrated in Fig. |[ We find 
that fl has to vary roughly like g^ in order that the spin 
gap is constant. Naively, one indeed expects that softer 
(i.e. with smaller Q) phonon modes are more effective to 
break the lattice symmetry. So, if one requires the spin 
gap to be constant, this effect has to be compensated by 
a smaller coupling g. 
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FIG. 8. Spin gap A/ J (in units of J) as a function of the 
magneto-elastic coupling g for a — 0.36, SI — 0.3J (o) and 
a — 0, ^ = 0.3 J (□). Horizontal marks indicate actual ex- 
perimental spin gap values. 

In Fig. ^ we observe that the coupling constant g(S7) 
is roughly 2.5 — 3 times smaller for CuGeOs than for 
NaV205 although the ratio of their spin gaps is only 1.5. 
This is an interesting consequence of the large frustration 
in CuGeOs. Indeed, a large a opens alone a (quite small) 
spin gap and, more importantly, amplifies the effect of the 
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spin-phonon coupling. This effect is even more drastic in 
the static model (||) where the dimerizations S = 0.014 
(CuGeOa) and 5 = 0.048 (NaV205) have a ratio of about 
4 lO. 
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FIG. 9. Frequency in units of J as a function of the 
magneto-elastic coupling g insuring a constant spin gap (see 
text) for a = 0.36 (o) and a = (□). 

The model (||) seems to describe accurately the spin- 
Peierls transition. Theoretical parameters have been de- 
duced from experiment and the ground state properties 
of the spin-Peierls phase have been established. We have 
provided evidences in favour of the dynamical breaking 
of the lattice periodicity with the simultaneous opening 
of the spin gap. Next, we shall study the dynamical prop- 
erties of this model. 



VI. DYNAMICAL PROPERTIES 

INS is a powerful experiment probing the momentum- 
dependence of the spin dynamics. INS has been per- 
formed on CuGeOa single crystals [ p^|jT^ and on NaV2 05 
powders |^ . It provides a direct measure of the dynami- 
cal spin-spin structure factor, 

S,,{q,cu) = ^ |(vI/„|5,((z)|vI/o)|25(cc> - i?„ + E„), (5) 



where jvl/o) is the (singlet) ground state of energy Eq and 
the sum is performed on all triplet excited states \'^n) (of 
energy En). Sz{q) is normahsed as i/VLJ2j G^p{i<lj)Sf ■ 
The INS spectrum can be easily computed by exact di- 
agonalization techniques Results on a 20 site chain 
are shown in Fig. |l^(a) for CuGeOa and in Fig. |l0|(b) 
for NaV205 with a frequency O = 0.3J. In both cases 
we observe a well defined q-dependent low energy struc- 
ture as for the static model (^. Its bandwidth (i.e. the 
energy at the maximum of the dispersion at g = tt/2) is 
typically tOmax ~ 1.1 J for CuGeOa and tOmax ~ 1.6 J 



for NaV205. This second value is very close to the 
DesCloizeaux-Pearson value of 7r/2 |38) of the Heisenberg 
chain in contrast to the case of CuGeOa which exhibits 
a large frustration. The ratio tOmax/J could therefore 
be considered as an additional accurate measure of the 
amount of frustration within the chain since the param- 
eter a alone determine approximately ujmax/ J ■ It is in- 
teresting to notice also that, in the case of a frustrated 
chain (CuGeOa), the upper limit of the continuum seems 
to be better defined. 

At low energy, the dimcrization gap leads to major 
differences with respect to the Heisenberg chain. First, 
there is no intensity for w < A. Secondly, the magnon 
branch is well separated from the continuum by a finite 
gap (see below) so that the magnon excitation can be 
interpreted as a spinon-spinon bound state p^ ]. This 
bound state was also found in the static model (^) for 
CuGeOs p3,pl and NaVaOg pTt]. 




FIG. 10. Szz{q,(jj) as a function of uj/ J calculated on a 20 
site chain for q = n^^, n = 0, • • ■ , 10; (a) CuGeOs parame- 
ters, a = 0.36, g = 0.109, Q. — 0.3J; (b) NaV205 parameters, 
a = 0, g — 0.270, fl — 0.3J. A broadening of the 5-functions 
e = 0.04 J was used. 

The dispersion relations of the magnon branch (o), the 
second excitation (□) and the upper limit of the con- 
tinuum (0) in the dynamical model are plotted in 
Fig. |ll|(a) for CuGeOa. The symbols correspond to ex- 
perimental results from Ref. [ p^ and filled symbols corre- 
spond to infinite size extrapolations at momenta q = it /2 
and g = vr. Similar dispersion relations are shown in 
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Fig. [TT|(b) for NaV205 and the position of the experimen- 
tal q = TT spin gap is indicated by an arrow. Note that 
we have expHcitly checked that the magnon branch is weU 
separated from the continuum. A finite size scahng anal- 
ysis of the energies of the two lowest triplet states (^M) 

(a-b) 
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is indeed possible at momentum q = 7r/2. Figs, 
clearly show that there is a finite gap between the first 
branch and the continuum as in the static model. It is 
consistent with the fact that the continuum corresponds 
to solitonic spin-1/2 excitations (or spinous) and that 
solitons and antisolitons can bind in pairs with momenta 
close to q = tt/2 Such a dougle gap feature was 

indeed observed experimentally p9| . 

It is important to notice that the dispersion relation is 
not symmetric with respect to g = 7r/2 in contrast to the 
case of a static dimerization. In fact, such a symmetry 
in the energy spectrum is due to the Bragg scattering 
resulting from the doubling of the unit cell. Since the 
dimerization appears only as a true phase transition in 
model (^), we expect that the symmetry of the spectrum 
with respect to tt/2 will only become exact in the thermo- 
dynamic limit. In the case of CuGeOa, our results are in 
very good agreement with INS experiments although fi- 
nite size effects are still important. In fact the agreement 
improves with increasing system size since the calculated 
magnon branch for q > it/2 shifts slightly to lower en- 
ergy when L grows (in order to be symmetric with the 
q < tt/2 part). Note also that energy scales are four times 
larger for NaV2 05 than for CuGeOa which could restrict 
INS experiments on NaV205 to low energy regions of the 
spectrum in the vicinity of g = tt. 

It is interesting to compare results for the spin dynam- 
ics obtained within the dynamical model to the ones ob- 
tained within the static model. Fig. |l2| shows the lowest 
triplet magnon branches and the next triplet excitations 
(in fact lower limits of the S ^ I continuum) for parame- 
ters suitable for CuGeOa. We do not explicitly show the 
comparison of the upper limits of the continua since the 
two curves obtained within the two models are almost 
indistinguishable. This is not surprising because higher 
energy excitations are only determined by the magnitude 
of the frustration and the coupling to the lattice plays a 
minor role here. At lower energy, the magnon branches 
of the two models look also very similar for q < tt/2 but 
some differences appear for q > tt/2 since, as explained 
before, the dispersion is not symmetric with respect to 
7r/2 in the dynamical model. This is simply due to larger 
finite size effects [0 occurring in model ^) related to the 
fact that the lattice periodicity is only spontaneously bro- 
ken. Once such finite size effects are taken into account 
we can safely conclude that the dispersions of the magnon 
branches of the two models in the thermodynamic limit 
are very close. Similarly, the discrepancies seen between 
the positions of the lower limits of the continua of triplet 
excitations are not relevant. Indeed, a detailed finite size 
scaling analysis at e.g. q = tt/2 reveals that the position 
of the two lower limits are in fact quite close (1.117 J for 
(^ to be compared to 1.118 J for (||). An exactly similar 
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FIG. 11. Momentum dependence of the first (o), second 
excitation (□) and upper limit of the continuum (0) on a 20 
site chain for the dynamical model (^. Filled symbols rep- 
resent extrapolations to infinite size (first (•) and second (■) 
excitations), (a) CuGeOa parameters (a — 0.36, g = 0.109, 
fl = 0.3 J). Experimental values (*) are taken from from 
Ref. Units on the right are in meV assuming J = WOK 
(13.8 meV). (b) NaV205 parameters {a = 0, g = 0.270, 
n = 0.3 J). Units on the right are in meV assuming J — 440K 
(37.9 meV). The arrow indicates the experimental value of the 
q — TV spin gap. 



comparison can be done for NaV2 05 (not shown). 
The spin static structure factor, 



Szziq) = J dujSzz{q,Lj) , 

which can be obtained in INS by integrating the spectrum 
over energy is plotted in Fig. ^ for CuGeOa (a = 0.36, 
g = 0.109, n = 0.3J) (o) and NaVzOg (•) (a = 0, g = 
0.270, O = 0.3J) for a 20 site chain. It is peaked near 
q = TT as a result of strong short range AF correlations. 
Indeed the width of the peak at q = tt is directly related 
to the inverse magnetic correlation length. Note however 
that Szz{tt) is slightly suppressed in CuGeOa compared 
to NaV205 because of the interchain frustration. In any 
case, the results are very similar to those obtained with 
the static dimerized model. The relative weights of the 
magnon peak in Szz{q, are also shown for CuGeOs (□) 
and NaV205 (■). Their behaviors versus q suggest that 
working in a range of momenta around q = O.Stt might be 
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FIG. 12. Momentum dependence of the two lowest triplet 
excitation energies in CuGeOs calculated on a 20 site chain for 
(i) the dynamical model (|) (a = 0.36, g = 0.109, D. = 0.3 J) 
(o) and (ii) the static model (|) (q = 0.36, S = 0.014) 
(□). Units on the right are meV assuming that J = 160A' 
(13.8 meV). 



more appropriate experimentally in order to have clearer 
evidences for the continuum. 



VII. CONCLUSIONS 



FIG. 13. Static factor structure for CuGeOa {a = 0.36, 
g = 0.109, D. = 0.3J) (o) and NaV205 {a = 0, g = 0.270, 
n = 0.3J) (•) calculated on a 20 site chain. Squares corre- 
spond to the weight of the lowest peak for CuGeOa (□) and 
for NaV205 (■). 



comments and IDRIS (Orsay) for allocation of CPU time 
on the C94 and C98 CRAY supercomputers. 

After completion of this wo rk, we learnt of a r elated 
work by A. W. Sandvik et al. ( :ond-mat /9706046 ) using 
a Quantum Monte Carlo approach. 



In order to describe one dimensional spin-Peierls com- 
pounds, a magneto-elastic (i.e. spin-phonon) coupling 
has been considered and is shown to be responsible for a 
dynamical and spontaneous breaking of the lattice peri- 
odicity followed simultaneously by the opening of a spin 
gap. The resulting symmetry-broken ground state is con- 
sistent with the existence of a frozen dimerization such 
as the one obtained in a mean-field treatment of the cou- 
pling to the lattice. We have used exact diagonalization 
techniques to calculate static and dynamical properties 
of this model. Controlled truncation procedures have 
been applied to the bosonic Hilbcrt space of the Hamil- 
tonian. By using a finite size scaling analysis, we have 
compared various physical quantities to the experimental 
ones (in the case of CuGeOa and NaV205) and we have 
determined a range of suitable parameters for the model. 
We predict that the spin-phonon coupling is 2 or 3 times 
larger in NaV2 05 than in CuGeOs. The INS spectrum 
calculated within this model is found to be qualitatively 
similar to the one obtained in the static model with a 
finite gap separating the magnon branch from the con- 
tinuum of triplet excitations above. 

D.A. acknowledges useful discussions with M. Albrecht 
and S. Capponi. We thank L. P. Regnault for communi- 
cating to us the data of Ref. J. Riera for valuable 
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